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Q ; Abstract 

p\J \ We have extended, from order 12 through order 25, the high-temperature series expansions (in 

P5 ' zero magnetic field) for the spin-spin correlations of the spin-S" Ising models on the square, simple- 

, ^ ■ cubic and body-centered-cubic lattices. On the basis of this large set of data, we confirm accurately 

-. . the validity of the scaling and universality hypotheses by resuming several tests which involve the 

^Sj ! correlation function, its moments and the exponential or the second-moment correlation-lengths. 
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I. INTRODUCTION AND CONCLUSIONS 

Moderate-length high-temperature (HT) expansions (through order 12) and low- 
temperature (LT) expansions for the spin-spin correlation function (sscf ) G{f, T; S) of the 
nearest-neighbor Ising models with general spin S were first computedi^ three decades ago 
on various lattices in 2D and in 3D. Motivations for the study of these models came not 
only from their direct phenomenological interest, but mainly from the conjecture^ that, in 
a given space dimension, the exponents characterizing the critical behavior are independent 
both of the lattice structure and of the spin magnitude S. This conjecture was the first 
step towards the modern notion of universality class. In the same years also the hypothesis 
of critical scaling^ was put forward. Many studie o^i^i^i^i'^i^i^i^°i^'^i^^i'^^ of the mentioned HT 
and LT series were devoted to test the validity and the main consequences of these basic 
hypothesea^ i^i^i^^i^^i^^i^^i^'^i^^i'^^ . Although the results sometimes were not as precise as was 
hoped, or covered only the 5=1/2 case, the scaling tests suggested that the critical sscf 
is a homogeneous function of appropriate variables, while the universality tests indicated 
that the critical indices and suitable combinations^^ of critical amplitudes are independent 
of the spin S and lattice structure. A few years later, the first substantial extension^iiS^ of 
HT Ising series in 3D (through order 21 on the body-centered-cubic(bcc) lattice only) did 
not make higher expansion coefficients available for the sscf, but only for its two lowest even 
moments and therefore various tests could not be repeated and updated. 

We are now resuming the HT part of those pioneering analyses in order to improve their 
extent and accuracy by taking advantage of our recent extensio n^^i^^i^^ from order 12 through 
order 25 of the HT expansions for the sscf of the Ising model with general spin S, in 2D 
on the square (sq) lattice and in 3D on the simple-cubic (sc) and the bcc lattices. From 
these data we have also derived series for related quantities, in particular for a variety of 
moments of the sscf, which are computed through order 25, and for the exponential (or 
'true') correlation length defined via the exponential decay of the sscf, which, however, can 
be extended only through order 19. For reasons of space we have not tabulated in this paper 
the series analysed, but have included them into our on-line library^^ of HT data for the 
spin-S" Ising model in order to make them more widely available for further study. Since this 
is the largest body of series data so far computed for these systems, we have already been 
studying other aspects of them in previous papers. In particular, in Ref .— we have accurately 
confirmed that the residual weak spin- dependence observed^- in lower order studies of the 
susceptibility exponent 7 and of the correlation-length exponent u in 3D on the bcc lattice, 
should not be ascribed to small violations of universality, but can be simply explained 
away as numerical inaccuracies due to expected non-negligible spin-dependent corrections 
to the leading scale behavior. Moreover we have tested the universality of several amplitude 
combinations obtaining similar results. In Ref.— an analogous survey of universal quantities 
was performed in 2D for the sq lattice case. Shorter series (but only for the 5 = 1/2 case) 
had been analysed in Refs.— . 

From the evidence presented here we can conclude that our HT data for the sscf have 
by now reached an extension sufficient to make the use of modern series-extrapolation tech- 
niques possible and generally reliable. Therefore we are able to exhibit more convincingly 
both in 2D and in 3D many expected properties related to scaling and universality also in 
some cases in which the old analyses led to inconclusive or not very precise results. 

The rest of the paper is organized as follows. In section II we shall outline the main 
features of the model, introduce our notations and conventions and very briefiy recall the 



scaling and universality properties expected for the sscf along with the corresponding tests 
discussed in full detail by the above cited papers^'^'SiSiiSiiiJ^. Therefore, in section III we 
can restrict ourselves to only a few comments on the numerical results. 

II. THE SPIN-S^ ISING MODELS 

The spin-S* Ising models with nearest-neighbor interaction are defined by the Hamiltonian: 

H{s} = -^ E ^(^1^(0 - ^E ^(^1 (1) 

where J is the exchange coupling, and s{f) = s^{f)/S with s^lf) a classical spin variable 
at the lattice site r, taking the 25" + 1 values —S, —S + 1, . . . , S — 1, S. The sum runs 
over all nearest-neighbor pairs of sites. For simplicity, the nearest-neighbor lattice spacing 
will be set equal to 1 everywhere. We shall consider expansions in the usual HT variable 
K = J/ksT where T is the temperature, fee the Boltzmann constant, and K will be called 
"inverse temperature" for brevity. In the critical region we shall also refer to the standard 
'reduced-temperature' variable t{S) = 1 — Tc{S)/T = 1 — K/Kc{S). 
We shall study the HT expansion of the (connected) sscf defined as 

G(f,T;5) = (s(0)s(r)), (2) 



In order to estimate numerically G{f, T; S*) as T — ^ Tc+ , we have allowed for its expected^*^^ 
behavior: in the 2D case 

G(r, T; S) ^ G{f, Tc] S) - E+{f; S)t{S)lnt{S) + ... (3) 

and in the 3D case 

G{r,T- S) ^ G{r,T,; S) - E+{r- S)t(5)i-" + . . . (4) 

Here E^{f;S) is the critical amplitude of the leading singular correction, a = 0.110(1)2^ 
denotes the critical exponent of the specific heat in 3D and the dots indicate higher order 
corrections. 

The correlation-function moment fin{T', S) of order n is defined as 

/x„(T;^) = EkT(s(0)^(^1)c (5) 

r 

(for n <0 the sum extends to r 7^ 0). 

The expected asymptotic behavior of finiT; 5*) as T — > Tc+ is 

UT; S) ^ m+(5)t(5)-(^+-)[l + a^iS)tiSf + ...]. (6) 

In 2D the exponent 6 of the leading singular correction is larger than unity, while in 3D 
a recent simultaneous study^^ of a set of models in the Ising universality class has suggested 
the very precise estimate 6 = 0.517(4). 

The scattering function, namely the Fourier transform of G{f, T; S) 

G{k,T-S)=J2exp[~tk-r\G{f,T;S) (7) 



for fc = yields the zero-field reduced susceptibility 

G(0, T; S) = /io(T; S) = x(T; S) = E(^(0)^(^1)c (8) 

The second-moment correlation-length is defined in d spatial dimensions by 

.2 /^ c-^ _ /^2(T;g) _ d\nG{k,T;S) 

^sm[^^^>- 2dx{T;S)~ dk^ ''=" ^^ 

For T > Tc the sscf is exponentially decreasing for large r and therefore following Ref.-^, 
beside the 'second-moment' correlation-length we can also define the inverse 'exponential' 
(or 'true') correlation-length in the direction e as 

kJT] S) = -limr^^ -ln\G(re,T;S)\ (10) 

r 

Since the singularity of G{k, T; S) closest to the real axis in the complex k plane is located 
at ±ik^{T; S), the exponential correlation-length can be obtained by solving recursively^ the 
eq. 

G{ikg,T;S)-' = 0. (11) 

Rather than working directly with kg{T; S) which is not an ordinary power series in K, it is 
expedient^ to form the quantity 



which is an ordinary power series in K. In Eq. ()12j) / is a geometrical factor depending on 
the unit vector e and on the lattice considered. In particular, if e is directed along a lattice 
axis, we have / = 1 for the sq and the sc lattices, while / = l/y/3 for the bcc lattice. 

So far, 3D data for this quantity were pubhshed exclusively for S = 1/2, and did not 
extend beyond order 15 in the sc lattice case^ or beyond order 10 in the bcc lattice case^. 
In 2D the HT expansion can be computed exactly^i^^ for S = 1/2, but no data have been 
published for S ^ \j1. In Ref.^^-, we have tabulated the expansion of ^|(T; 5*) through order 
19 for e directed along a lattice axis, in the case of the sq, sc and bcc lattices and with 
S= 1/2,1,3/2,2,5/2,3,00. 

In order to avoid possible confusion, it should be pointed out that in Ref.^ our ^ was 
denoted by A2(e), while the symbol ^^ was used to denote fcgT^. Our notation might be 
more suggestive since our ^| compares very closely with ^^^. Indeed, the true and second- 
moment correlation-lengths are almost identical in magnitude above the critical temperature. 
In particular on the sq lattice, when e is directed along a lattice axis, the HT expansion 
coefficients of and Q.^ coincide through sixth order for 5* = 1/2 , through fourth order 
for S = \ and through second order for higher values of the spin. In 3D, in the sc lattice 
case, the expansion coefficients of and Q.^ coincide through seventh order for S = 1/2, 
through fifth order for S = 1 and through third order for higher values of S. In the case of 
the bcc lattice, the expansion coefficients coincide through third order for all values of the 
spin. Moreover, up to the maximum order of our computation, the noncoinciding coefficients 
differ by less than 0.1%. 



The two correlation lengths ^g and ^sm are expected to share the same critical exponent 
V so that their asymptotic behavior when T — > Tc+ can be written as 

is^iT- S) ^ ft^{S)t{Sr^[l + at^{S)t{Sf + ...] (13) 

and 

UT: S) ^ r (5)t(5)-^[l + a^{S)t{SY + ..] (14) 

Here the critical amplitude f'^{S) is independent of e, since the sscf becomes spherically 
symmetric near the critical point. The ratio 

QtiS) = /-*-(5)//l(^) (15) 

is a universal combination of critical amplitudes^S i.e. it is expected to depend only on the 
lattice dimensionality d, but not on the spin S or the lattice structure. 

For T ^ Tc + and in zero magnetic field, the sscf is expected to exhibit the asymptotic 

stmcijurG 

Gir,T; S) ^ {l/ry-'+^MS)Do{Ci{S)r/^,UT; S)) + . . . (16) 

when both r and C,sm are much larger than the lattice spacing (with arbitrary r/^sm)- Eq. p6p 
together with these assumptions on r and ^sm, is usually referred to as the "strong-scaling 
hypothesis" (while it is called the "weak-scaling hypothesis" , if its validity is restricted to 
the r ^ oo limit with fixed r/^gm)- In Eq. (jl6|) . r] is the critical exponent describing the 
decay of the sscf at the critical point, Do(a;) is called the critical scaling function, ^4^(5*) and 
Ci{S) are scale factors. The dots indicate subcritical corrections proportional to a positive 
power of some irrelevant field. The scaling function Do(x) is expected to be universal: its 
structure does not depend on the particular model under study provided that it belongs to 
a given universality class. On the contrary the scale factors ^^(5'), Ci{S) depend on the spin 
and the lattice /. The validity of the asymptotic structure Eq. ()16|) was verified analytically^ 
for the spin 5=1/2 Ising model in 2D. 

For the scattering function G{k,T]S) the analogous scaling form as T ^ Tc+ can be 
written as 

Gik,T;S)^A[iS)t{S)-^I)',{c[iS)eesJ)+--- (17) 

If the scale factors A[{S),Ci{S) are specified adopting the normalization conditions 



d;(0) = 1; (^).^o = -1- (18) 



one can write^ as fc -^ 



G(0, T; S)/Gik, T; S) = l/g^k^smiT; S)) = (19) 

1 + eL(T; S)k' - S4(T, S)CmiT; S)k' + S6(T, 5)eL(T; S)k' + 0{k^) 

where the function g+{k^smiT] S)) is universal and thus the quantities 

S,(T,5) = c /^^^',^^^°(f^^ -l (20) 

Se(T, S) = e /'^^^'^y^^'^) _ 2^,iT, S) - 1 (21) 

^''21-' 1 ^) 



with C4 = 1/4 and cg = 1/36 for d = 2, while C4 = 3/10 and ce = 3/70 for d = 3, have finite 
universal values as T — >■ Tc+. 

All T,2n(T, S), as well as the difference Q — ^g„, vanishS. in the mean-field related approx- 
imations. Therefore the magnitudes of these quantities at the critical point can be taken as 
a measure of the deviation of a given system from gaussian behavior, which turns out to be 
very small on the HT side of the critical point. 

More generally, it was observed^ that the scaling hypothesis Eq. (jl6p implies that, at the 
critical point, the ratios 

p (rp- Q^ _ l^'nriyT] S)iXn{T; S) , . 

IJ,r{T; S}iJ,s{T; S) 

with m + n = r + s, are universal. These ratios are dominated by the critical singularity 
also for negative values of the indices m, n, r, s provided that each index exceeds —2 + 77, as 
follows from Eq.®. 

Finally, the determination of the amplitude E~^{f; S) of the leading singularity of the 
sscf (see Eqs.dS)) and (JD)) gives another opportunity to perform universality and scaling 
tests. In order that the structure of Eqs.Q and (@]) be compatible^^ with the strong-scaling 
hypothesis Eq. ()16|) . the amplitudes E^{f; S) must scale as r'' with C = (1 — a)/h' + 2 — d — ri, 
namely 

E+{r; S) ^ E+{Sy (23) 

for large enough r, independently of the spin and the lattice structure. In 2D the value 
( = 0.75 is expected, while in 3D , adopting our recent estimate^^ of the values of the 
correlation length exponent u = 0.6299(2) and of the exponent rj = 0.036(1), we should have 
( = 0.3765(10). 



III. NUMERICAL RESULTS 

Let us first observe that, due to the leading singular corrections in Eqs.© and (JH), 
whose amplitudes E~^ (r ; S) grow with r as indicated by Eq. (PH|) , determining accurately 
G{f,Tc; S) ( as well as E~^{f; S) itself ) is a rather delicate matter for which it is crucial to 
rely on sufficiently many expansion coefficients. We should also consider that the number 
of non-trivial coefficients in our series decreases with increasing r, and correspondingly the 
precision of our estimates of G{f,Tc] S) (and E^{r, S)) deteriorates. In Refs .^1^^1^^ , due to 
the small number of coefficients available at that time, a generalized Neville extrapolation 
of the partial sums had to be used for determining G{f,T] S) and E~^{f; S) in the vicinity 
of Tc. Taking advantage of our new series, we can now improve substantially the numerical 
resummation of the HT series by resorting to first- or second-order inhomogeneous differen- 
tial approximant32i(DA's) biased with Kc{S). (Here and in what follows we have adopted 
the values of the critical temperatures tabulated in Refs .^^i^^ .) It does not come as a sur- 
prise that our procedures are slightly less efficient in 2D than in 3D, probably due to the 
presence of logarithms in the leading correction terms to the critical asymptotic behavior 
Eq.Q and also that in 3D the bcc lattice series always yield the most accurate results. If 
we restrict to 1. < r < 6. the relative uncertainty of our estimates of the critical sscf should 
generally remain well below 1%. In 2D this can be guessed by comparing the estimates of 
G{f,Tc; S) obtained from our series 0(Ar^^) with the known exact results in the sq lattice 
case^»^ for S' = 1/2 and safely assuming that the precision does not deteriorate too fastly 
when higher values of S are considered. In 3D no exact results are available, but the HT 

6 



series for the nearest-neighbor correlation function was recently extended^ through order 
45 in the sc lattice case for S = 1/2. Therefore, in this case, we are able to compare our 
estimate at order 25 with the result obtained by applying the same numerical procedures 
to the series extended through order 45. (It would be very interesting if the improved 
finite-lattice technique devised for this remarkable calculation could be generalized as effec- 
tively beyond first-neighbor correlations and to general S.) We should also mention that a 
completely consistent alternative estimate of the critical sc-lattice nearest-neighbor sscf has 
been obtained^^ in a recent high-precision MonteCarlo study. For other values of r in the sc 
lattice case and in the bcc lattice case our results can only be compared with calculations— 
using the old series 0{K^^). Table I lists our estimates of G{f,Tc] S) with their apparent 
uncertainties for a small sample of values of f and S. Previous estimates of the critical sscf 
from shorter series, which are available only for S = 1/2, are shown for comparison in the 
first column, labelled [5* = 1/2], of this table. In Figs lll2l51 we have plotted our estimates of 
lnfG(r, Tc, S)) vs ln(r) for 1 < r < 5 with S = 1/2, 1, 3/2, 2 in the cases of the sq, sc and bcc 
lattices respectively. We have also shown by continuous lines the results of one-parameter 



TABLE I: Our estimates of the critical-point values G{f,Tc; S) of the spin-spin correlation function 
for the nearest-neighbor Ising models with spin S = 1/2, 1, 3/2, 2, cxd on the sq, sc and bcc lattices. 
For comparison with our results, the first column labelled [S = 1/2] shows the only available 
previous estimates. In the case of the sq lattice, we have cited the exact^i^^ values. In the case of 
the nearest-neighbor correlation (r = (1, 0, 0)) on the sc lattice, we have reported in the first column 
the estimate obtained using our numerical procedure with the series 0{K^^) of Refi^ and have 
also cited a valued obtained in a recent high-precision MonteCarlo simulation. In the remaining 
cases we have reported the estimates of Refi^ obtained from series 0{K^'^). We are not aware of 
previous calculations for S > 1/2. 



Lattice f 


[S=l/2] 


S=l/2 


S=l 


S=3/2 


S=2 


S=cx) 


sq (1,0) 


0.707107..'^ 


0.7071(1) 


0.5806(3) 


0.517(1) 


0.481(1) 


0.338(1) 


(1,1) 


0.636620..'^ 


0.6366(1) 


0.5207(4) 


0.463(1) 


0.431(1) 


0.303(1) 


(2,0) 


0.594715.." 


0.5947(2) 


0.486(1) 


0.433(1) 


0.402(1) 


0.282(2) 


(2,1) 


0.573159.." 


0.573(1) 


0.467(1) 


0.417(2) 


0.387(2) 


0.272(2) 


(2,2) 


0.540380.." 


0.540(1) 


0.442(1) 


0.393(2) 


0.365(2) 


0.256(4) 


sc (i,o,o: 


0.330200(5)^ 


0.33020(6) 


0.24203(6) 


0.20756(6) 


0.18918(6) 


0.12886(6) 


(i,o,o; 


0.33017(3)^= 












(i,i,o: 


0.208(2)"' 


0.2086(1) 


0.1529(1) 


0.1311(1) 


0.1194(1) 


0.08141(5) 


(i,i,i: 


0.164(4)"' 


0.1633(1) 


0.1197(1) 


0.1027(1) 


0.0936(1) 


0.0638(1) 


(2,o,o: 


0.162(4)"' 


0.1608(2) 


0.1178(2) 


0.1010(2) 


0.0921(2) 


0.0627(1) 


(3,o,o: 


0.104(7)'^ 


0.1017(3) 


0.0746(2) 


0.0639(2) 


0.0581(2) 


0.0396(1) 


bcc (i,i,i: 


0.2735(7)^^ 


0.27265(5) 


0.19653(5) 


0.16763(5) 


0.15243(5) 


0.10341(5) 


(2,o,o: 


0.200(2)-^ 


0.19971(5) 


0.14394(5) 


0.12278(5) 


0.11165(5) 


0.07575(5) 


(2,2,0: 


0.157(2)-^ 


0.15627(5) 


0.11269(5) 


0.09614(5) 


0.08743(5) 


0.05934(5) 


(3,i,i: 


0.129(3)'^ 


0.12751(5) 


0.09193(5) 


0.07843(5) 


0.07132(5) 


0.04839(5) 


(2,2,2: 


0.131(3)'^ 


0.12914(5) 


0.09315(5) 


0.07945(5) 


0.07224(5) 


0.04903(5) 




"Reference^ 


^Reference^ 


'^References '^ReferenceS 





fits to the leading asymptotic behaviors lnfG(r, Tc; 5*) j ~ c{S) — {d — 2 + ?7)ln(r) expected 
for large r. We have taken only c{S) as a free parameter and fixed rj = 0.25 in 2D and 
T] = 0.036 in 3D. 

Both in 2D and in 3D, we have estimated also E^{f] S) from the amplitude of the singu- 
larity of the second temperature derivative of G{f,T] S), again using inhomogeneous first- 
and second-order DA's biased with Kc{S) and a. Our estimates of E'^{f]S) for a small 
sample of values of r and S are shown in Table II. They are compared with the exactly 
known values^ for S = 1/2, in the case of the sq lattice, or with a few old estimates-^ from 
shorter series, in the case of the sc and bcc lattices. A comparison with the exact results in 
2D and with our estimate using the mentioned high-order calculation in the sc lattice^^ for 
S = 1/2, still suggests that, for all values of S, the relative accuracy of our estimates should 
not be generally worse than 1%. 



TABLE II: Amplitudes E^{f;S) of the leading singular correction of the sscf near the critical 
point for the nearest-neighbor Ising models with spin S = 1/2, 1, 3/2,2, oo on the sq, sc and bcc 
lattices. For comparison with our results, the first column of the table labelled [S = 1/2], shows 
the available estimates from other sources. In the case of the sq lattice, the exact values are taken 
from Refi^. In the case of the nearest-neighbor correlation on the sc lattice (r = (1,0,0)), we 
have reported in the first column our estimate obtained from the series 0{K'^^) of Ref.^^. In the 
remaining cases, whenever available, we have quoted the estimates of Ref.^^ obtained from series 
0{K^'^). We are not aware of other published calculations for S > 1/2. 



Lattice 



[S=l/2] 



S=l/2 



S=l 



S=3/2 



S=2 



=CXD 



sq 



sc 



bcc 



(1,0) 


0.561100.." 


0.562(1) 


0.621(1) 


0.623(1) 


0.613(1) 


0.484(1) 


(1,1) 


0.793515.." 


0.794(1) 


0.819(1) 


0.812(1) 


0.794(1) 


0.616(1) 


(2,0) 


1.0103348.." 


1.01(1) 


1.02(1) 


1.01(1) 


0.987(2) 


0.759(2) 


(2,1) 


1.120022.." 


1.11(1) 


1.13(1) 


1.11(1) 


1.08(1) 


0.826(2) 


(1,0,0) 


2.252(5)^^ 


2.27(2) 


2.16(2) 


2.03(2) 


1.93(2) 


1.42(2) 


(1,1,0) 


2.38(2)^= 


3.01(2) 


2.72(2) 


2.52(2) 


2.38(2) 


1.72(2) 


(1,1,1) 


2.86(4)'= 


3.40(2) 


3.03(2) 


2.78(2) 


2.62(2) 


1.90(2) 


(2,0,0) 


3.16(6)^ 


3.53(2) 


3.14(2) 


2.88(2) 


2.71(2) 


1.95(2) 


(3,0,0) 




4.36(2) 


3.79(2) 


3.45(2) 


3.24(2) 


2.33(2) 


(1,1,1) 


2.010^= 


2.325(5) 


2.167(5) 


2.022(5) 


1.917(5) 


1.401(5) 


(2,0,0) 




2.707(6) 


2.442(6) 


2.256(6) 


2.129(6) 


1.545(6) 


(2,2,0) 




3.126(6) 


2.767(6) 


2.535(6) 


2.384(6) 


1.720(6) 


(3,1,1) 




3.41(1) 


2.98(1) 


2.72(1) 


2.55(1) 


1.83(1) 


(2,2,2) 




3.44(1) 


3.01(1) 


2.74(1) 


2.57(1) 


1.84(1) 



"Reference^ ^Reference^^ '^Referenceii 



In Figs. I4I5I6I for S = 1/2,1,3/2,2 we have plotted, \n(E^{f;S)) vs ln(r) in the case 
of the sq, sc, and bcc lattices respectively. For r > 4.5 in the case of the sc lattice and 
r > 6. in the case of the bcc lattice, we have not reported any estimates of E~^{'r;S), 
because the available nontrivial HT coefficients of the sscf are not sufficiently many to allow 
estimates at the level of precision above mentioned. In these figures we have also represented 
by continuous lines the results of one-parameter fits to the leading asymptotic behaviors 
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Inf £'+(r; 5") j ~ b{S) + C^n{r) expected for large r. We have taken for ( the expected values 
( = 0.75 in the 2D case and ( = 0.3765 in the 3D cases, while the free parameters b{S) 
have been determined using in the fits only the data with r ^ 1.8. Indeed, our new data 
show visible deviations from asymptotic scaling for sufficiently small r, particularly so in 
the case of the sc lattice, but the asymptotic consistency with the strong-scaling hypothesis 
Eq. fl23p is good. The behavior of E^{f; 1/2) as a function of r was first studied in Ref- 
using series 0{K^^) for the face-centered cubic lattice. In that analysis both ( and b{l/2) 
were determined by a two-parameter fit of the numerical results to the leading asymptotic 
behavior under the very simple assumption that the corrections to scaling are negligible 
even for 'not very large' r. (As we have indicated above, our new data show that such 
a strong assumption is untenable.) The authors of Ref.— concluded that ( = 0.47(6), an 
estimate in sharp disagreement with the value ( = 0.33(1) expected from the exponent 
values z/ = 0.638(2) and t] = 0.041(6) generally accepted at that time. A few years later, 
for S = 1/2, the somewhat lower estimate ( = 0.39(4) was obtaineciii from an analysis of 
LT expansions in powers of m = exp{—2K) up to order u^^ and u^^, on the sc and the bcc 
lattices respectively. In this latter study, however, a third fit parameter had been introduced 
in order to allow for small corrections to the asymptotic scaling behavior of E~^ {-r; S) . Also 
this estimate of ( did not agree with the value expected at that time, but is quite compatible 
with the presently preferred value. 

Before any strong confidence in the results of such two- or three-parameter fits can be 
justified, we believe, however, that the HT series should be further extended in order to 
enlarge significantly the range of values of r for which E~^{f]S) can be determined with 
sufficient accuracy. 

Having tabulated a wide sample of estimates of G{f,Tc] S) and E^{f] S) with some im- 
provement both in the extent and the accuracy, with respect to the very few estimates 
available in the literature, we are now in the position to exhibit more directly the scaling 
property by examining the near-critical sscf in the r-space. For T ^ T^ -|- 0, as suggested by 
Eq. (fTH|) . by a proper choice of the scale factors Ai{S) and Ci{S), we should be able to plot 
the quantities 

r'i-^+r,G{f,T; S) ^ A{S)Do(Q{S)r/^,UT; S)) (24) 



vs. r/^sm in such a way that the curves, associated to various values of S and to different lat- 
tices, collapse on each other. In Fig|7l we have plotted \n(r'^^'^^^'G{f,T; S)) vs. r/^sm{T] S) 
in the case of the sq lattice taking r = (2, 0) and S = 1/2, 1, 3/2, 2. Our data points refer to 
the range of temperatures for which 1.5 ^ ^sm ~ 200. Fig. |H1 shows the analogous plot for 
ln(r'^~^+''G(r, T; S)) vs. r/^sm in the case of the sc and bcc lattices. Here we have taken 
r = (4, 0, 0) and S = 1/2, 1, 3/2, 2 and have plotted data in the range of temperatures for 
which 2.7 ^ ^sm ~ 400. Completely consistent results are obtained also for other choices of 
r. As already observed, within these limitations, the present length of the HT series appears 
sufficient to obtain reliable estimates and our results are consistent with the strong-scaling 
hypothesis to a good approximation. 

The very small mismatch of the curves in the extreme regions r/E,siri ^ 1 or r/^sm 3> 
1 which can still be observed is related: i) to the fact that the scaling property has an 
asymptotic character, while in practice the size of r still cannot exceed a few lattice spacings 
if we want to use a decent number of expansion coefficients in the estimate of G{f, T; S), ii) 
to the residual infiuence of the subcritical corrections. 

We can further test the universality properties of the sscf in the /c-space, namely the 
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critical scattering function, by simply showing that T^^iTc, S) and T,q(Tc; S) are independent 
of S and of the lattice structure. Also these quantities are calculated by first- and second- 
order DA's biased with Kc{S). Since higher-order moments of the sscf (in which the less 
accurately known correlations between distant spins are weighted much more than those 
between near spins) enter into the definitions eq. (j2m) and eq. (j2ip . the convergence of the 
extrapolations is not expected to be very fast, particularly so in the cases of the sq and sc 
lattices. We should also consider that H^lTc, S) is the very small difference between unity 
and the critical value of some multiple of a ratio of moments of the sscf, so that a very high 
accuracy in the estimate of the latter is needed to achieve even a relatively modest precision 
for S4(Tc; S). The same remark applies also in the case of Tiq{Tc] S). In Table UTTl we have 
collected our estimates of E4(Tc; S) and T,q{Tc] S) in the case of the sc, sq and bcc lattices for 
5* = 1/2, 1, 3/2, 2, cxD. We have also reported a few previous estimates'^'^'^^ from the existing 
literature. 

In the case of the sq lattice our data suggest the final estimates T^^^Tc, S) = 7.8(3) x 10~^ 



TABLE III: Estimates of J^iiTc, S) (see Eq.d^n))) and T.q{Tc; S) (see Eq.dHI) in the case of the sq, 
sc and bcc lattices for various values of S. For comparison, we have also reported a few previous 
estimates listed in Ref i^'^i^^ , indicating the method of calculation and the expected uncertainty, 
when available. Only for convenience, we have listed in the S = 1/2 column also the results of 
the renormalization group and the optimized continuous spin calculations which, of course, do not 
refer to spin S = 1/2. 



Quantity 




Lattice 


S=l/2 


S=l 


S=3/2 


S=2 


S=cx) 


^T, 


S) X 10^ 




sq 


7.8(3) 


7.9(3) 


7.9(3) 


7.6(3) 


7.5(3) 


U{T, 


S) X lO^(Exact)'^ 


sq 


7.936796. . . 










^6{Tc 


S) X 10^ 




sq 


1.1(1) 


1.1(1) 


1.0(1) 


1.1(1) 


1.0(1) 


^6{Tc 


S) X 105(Exact)'^ 


sq 


1.095991... 










S4(Te 


S) X 10^ 




sc 


3.76(8) 


3.9(2) 


3.77(8) 


3.75(8) 


3.7(2) 


^4{Tc 


S) X 10^ 




bcc 


3.75(5) 


3.74(5) 


3.76(5) 


3.76(5) 


3.77(5) 


S6(T, 


S) X 10^ 




sc 


1.0(2) 


.9(2) 


0.9(2) 


0.8(2) 


0.7(2) 


S6(T, 


S) X 10^ 




bcc 


0.9(1) 


0.86(5) 


0.85(5) 


0.85(5) 


.85(5) 


S4 X lO^pT]^ 




sc 


3.0(2) 










S4 X 10^ [HT]^ 




sc 


5.5(15) 










S4 X 10^ [HT]^ 




bcc 


7.1(15) 










Se X lO^pT]'' 




sc 


0.5(2) 










Se X 105[HT]'= 




sc 


0.5(2) 










Se X 105[HT]'= 




bcc 


0.9(3) 










S4 X 10^ [opt.cont. 


spin]^ 


sc 


3.90(6) 










Se X 10^ [opt.cont. 


spin]^ 


sc 


.88(1) 










S4 X 10^ [e-expans 


f 




3.3(2) 










Sg X 10^ [e-expans 


? 




0.7 










S4 X 10^ [g-expans 


? 




4.0(5) 










Sg X 10^ [g-expans 


]' 




1.3(3) 











" Reference2£ ^Reference^y^ '^Reference^ 
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and I1q{Tc; S) = 1-1(1) x 10 ^, independently of S and in reasonable agreement with the 
high-precision determinations^^ of J^^iTc, 1/2) = 7.936796. . . x 10"^ and of ^eiTc, 1/2) = 
1.095991 ... X 10~^ obtained by numerical integration of the analytically known3^ sscf of 
the S = 1/2 model in 2D. In 3D our results for the bcc lattice show a definitely smaller 
uncertainty than for the sc lattice. They suggest the final estimates S4(Tc; S) = 3.8(1) x 10~^ 
and T,q(Tc; S) = 0.9(1) x 10~^ independently of the spin 5" and lattice structure. Our results 
are therefore consistent with the corresponding estimates in the literature, in particular with 
the values E^iTc) = 3.90(6) x 10"^ and S6(Tc) = 0.88(1) x 10"^ obtained^^^ optimizing the 
parameters of a continuous-spin model, under the assumption of universality. Let us also 
mention that renormalization group calculations^i in the e-expansion scheme to third order 
yielded the estimates S4 = 3.3(2) x 10~^ and Eg = 0.7 x 10~^, while, in the coupling-constant 
expansion scheme to fourth order, the corresponding results were S4 = 4.0(5) x 10~^ and 
Eg = 1.3(3) X 10"^ 

TABLE IV: Estimates of the moment ratios Rm,n;r,s{Tc; S) (see Ea. (|22p l in the case of the sq, sc 
and bcc lattices for various values of S. 



i^Tn,n;r,s 




Lattice 


S=l/2 


S=l 


S=3/2 


S=2 


S=oo 


-^0,l;l/2,l/2 




sq 


1.1641(1) 


1.1642(1) 


1.1642(1) 


1.1642(1) 


1.1641(1) 


-^0,l;l/4,3/4 




sq 


1.1211(1) 


1.1211(1) 


1.1211(1) 


1.1211(1) 


1.1210(1) 


-^-3/4,l/4;-l/4,- 


-1/4 


sq 


1.299(1) 


1.300(1) 


1.301(1) 


1.300(1) 


1.301(1) 


-^-l,-l/2;-3/4,- 


-3/4 


sq 


1.121(5) 


1.124(4) 


1.124(4) 


1.125(4) 


1.126(4) 


-Ro,l;l/2,l/2 




sc 


1.1320(1) 


1.1320(2) 


1.1319(1) 


1.1319(2) 


1.1319(2) 


-^0,l;l/2,l/2 




bcc 


1.1320(1) 


1.1319(1) 


1.1319(1) 


1.1319(1) 


1.1319(1) 


-^0,l;l/4,3/4 




sc 


1.0977(2) 


1.0977(2) 


1.0976(2) 


1.0976(2) 


1.0976(2) 


-^0,l;l/4,3/4 




bcc 


1.0977(1) 


1.0976(1) 


1.0976(1) 


1.0976(1) 


1.0976(1) 


-Rl/2,l/2;l/4,3/4 




sc 


0.9697(2) 


0.9697(2) 


0.9698(2) 


0.9697(2) 


0.9697(2) 


-^l/2,l/2;l/4,3/4 




bcc 


0.9697(1) 


0.9697(1) 


0.9697(1) 


0.9697(1) 


0.9698(1) 


-R-l,-l/2;-3/4,- 


-3/4 


sc 


1.084(1) 


1.084(1) 


1.084(1) 


1.084(1) 


1.083(1) 


-R-l,-l/2;-3/4,- 


-3/4 


bcc 


1.083(1) 


1.083(1) 


1.083(1) 


1.083(1) 


1.083(1) 



The results of our analysis of the universal ratios Rm,n;r,s(Tc; S) are reported in table 
IIVI They also show independence of the spin and of the lattice structure within a good 
precision. Our series-extrapolation procedure, based on first- and second-order DA's uses 
only our estimates of Kc{S) and does not need to be biased also with 7 and z/ as it was 
necessary in the generalized Neville procedure^J^ employed with the short series of Ref.— . 
Considering that the values 7 = 1.25 and z/ = 0.625 (or z/ = 0.638) of the exponents accepted 
at the time of that study are somewhat different from the currently preferred ones and that 
the extrapolations are very sensitive to those values, a comparison with the numerical results 
of Ref.— has little meaning. 

Finally, we have tested both in 2D and in 3D the spin independence of the ratio Q^{S) 
defined by EqUHl In 3D also the lattice independence of Q^{S) can be tested. 

In 2D, on the sq lattice, the non-trivial expansion coefficients of the ratio 
^giT; S)/^'^^(T; S) are not sufficiently many and their behavior is not smooth enough to yield 
very accurate results. Therefore our best estimate of Qt{S) (by first-order DA's biased with 
Kc{S)) cannot be more precise than Q^{S) = 1.0004(2), independently of S. Our rough es- 
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timate is, however, consistent with the more accurate determination (5/(1/2) = 1.000402 . . . 
obtained in the S = 1/2 case in which, as aheady indicated above, very long series are 
available^i for ^^^(r; 1/2), while ^j{T- 1/2) is exactly known^-^s^ 

In 3D we can use both first- and second-order DA's biased with Kc{S). The very smooth 
bcc lattice series yield the most accurate results. Our final estimate is Q'^iS) = 1.000200(3), 
independently of 5* and of the lattice structure. So far, this ratio could be computed'^^ only 
for S = 1/2 from a 15 term series on the sc lattice, with the result Qt{l/2) = 1.000125(50). 
A more precise estimate^ Qt = 1.000199(3) was obtained indirectly (and assuming uni- 
versality), from optimized HT series for a continuous-spin model on the sc lattice. Within 
the renormalization group approach^^, the estimate Qt = 1.000160(20) was obtained in the 
e-expansion to third order, while the coupling-constant expansion technique to fourth order 
gave Q^ = 1.000205(30). 
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FIG. 1: Fig.l. Estimates of G{f, Tc] S) on the sq lattice. The meaning of the symbols is as follows. 
Triangles: S = 1/2, squares: S* = 1, rhombs: S = 3/2, circles: S* = 2. The spin S points are 
shifted vertically by the quantity 1/2 — S* in order to make the figure more legible. The continuous 
lines represent fits to the leading asymptotic behaviors lnG(r, Tc;S') ~ c{S) — 77ln(r) expected for 
large r. We have taken c{S) as a fit parameter and fixed rj = 0.25. 
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FIG. 2: Fig. 2. Estimates of G{r,Tc]S) on the sc lattice. The meaning of the symbols is the 
same as in Fig^The spin S points are shifted vertically by the quantity 1/2 — 5 in order to make 
the figure more legible. The continuous lines represent fits to the leading asymptotic behaviors 
lnG(f, Tc, S) « c{S) — (1 + r/)ln(r) expected for large r. We have taken c{S) as a fit parameter and 
fixed T] = 0.036. 



15 



CD 
C 




LO _ 



0.0 



0.5 



1.0 



1.5 



FIG. 3: Fig. 3. Estimates of G{f, T^, S) on the bcc lattice. The meaning of the symbols is the 
same as in Fig^ The spin S points are shifted vertically by the quantity 1/2 — 5 in order to make 
the figure more legible. The continuous lines represent fits to the leading asymptotic behaviors 
lnG(f, Tc, S) « c{S) — (1 + r/)ln(r) expected for large r. We have taken c{S) as a fit parameter and 
fixed T] = 0.036. 
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FIG. 4: Fig. 4. Estimates of E~^{r\ S) on the sq lattice. The meaning of the symbols is the same 
as in FigQ The spin S points are shifted vertically by the quantity 1/2 — S* in order to make 
the figure more legible. The continuous lines represent fits to the leading asymptotic behaviors 
ln£^+(f; S) « b{S) + Cln(r) expected for large r. We have taken 6(5) as a fit parameter and fixed 
C = 0.75. 
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FIG. 5: Fig. 5. Estimates of E^{f;S) on the sc lattice. The meaning of the symbols is the same 
as in FigQ The spin S points are shifted vertically by the quantity 1/2 — S* in order to make 
the figure more legible. The continuous lines represent fits to the leading asymptotic behaviors 
ln£^+(f; S) « b{S) + Cln(r) expected for large r. We have taken 6(5) as a fit parameter and fixed 
C = 0.3765. 
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FIG. 6: Fig. 6. Estimates of E^{f; S) on the bcc lattice. The meaning of the symbols is the same 
as in FigQ The spin S points are shifted vertically by the quantity 1/2 — S* in order to make 
the figure more legible. The continuous lines represent fits to the leading asymptotic behaviors 
ln£^+(f; S) « b{S) + Cln(r) expected for large r. We have taken 6(5) as a fit parameter and fixed 
C = 0.3765. 
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FIG. 7: Fig. 7. The logarithm of the scaUng function r '^~^'^G{f,T; S) vs. x = r/^sm{T; S) in the 
case of the sq lattice. The data represent the ssfc's with f = (2,0) and S = 1/2,1,3/2,2 in the 
range of temperatures for which 1.5 < ^(T; S) ^ 200. The meaning of the symbols is the same as 
in Fig^ 
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FIG. 8: Fig. 8. The logarithm of the scaUng function r'^~'^^^G{r,T; S) vs. x = r/^sm{T; S) in the 
case of the sc and the bcc lattices. For both lattices the data represent the ssfc's with r = (4, 0, 0) 
and S = 1/2, 1, 3/2, 2 in the range of temperatures for which 2.7 ^ (.{T; S) ^ 400. The meaning of 
the symbols is the same as in Fig^for the sc lattice case. For the bcc lattice data we have used 
full symbols of the same shape. 
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